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EXECUTIVE  SUMMARY 

This  study  of  corrections  to  propagation  predictions  based  on 
ordinary  ray  theory  is  motivated  by  the  problem  of  low-frequency 
ambient  noise  caused  by  surface  shipping. 

Shipping  noise  is  caused  by  many  sound  sources,  widely  dis- 
tributed in  range  and  azimuth  from  any  receiver.  The  sources  - 
propellers  and  machinery  — are  near  the  sea  surface,  at  depths 
where  large  gradients  of  the  sound  speed  are  often  found.  These 
circumstances  can  lead  to  significant  errors  in  calculations  of 
low-frequency  sound  propagation  based  on  ordinary  ray  theory. 

The  objective  of  this  study  is  to  quantify  the  significant 
errors  in  ray  theory  calculations,  and  to  indicate  when  they  arise 
and  how  corrections  may  be  made  in  the  calculation  of  low-frequency 
shipping  noise.. 

The  most  significant  error  concerns  the  acoustic  source 
strength.  When  the  source  ship  is  in  a region  having  large  nega- 
tive values  of  sound-speed  gradient,  ordinary  ray  theory  under- 
estimates the  sound  power  injected  into  ray  bundles  with  shallow 
D/E  angles.  Errors  of  6 dB  and  more  are  found  in  realistic 
circumstances . 

Correction  for  this  error  is  readily  achieved  by  using  an 
effective  source  strength  that  depends  on  frequency,  D/E  angle, 
and  the  sound  speed  profile  in  the  surface  layers.  The  effective 
source  strength  amounts  to  a wave  theoretical  modification  to 
Lloyd's  mirror  formula,  which  has  often  been  used  in  ray  theory 
to  describe  the  reduction  in  source  strength  due  to  a source 
being  located  near  the  ocean  surface. 


A 


Report  No.  3326  Bolt  Beranek  and  Newman  Inc. 

A second  error  concerns  the  path  of  rays.  Ordinary  ray  theory 
overestimates  the  range  interval  for  which  shallow  rays  dwell  in 
a surface  layer  having  a large  sound-speed  gradient . The  error 
occurs  on  every  cycle  to  the  surface,  and  is  cumulative  with  range. 
In  realistic  circumstances,  errors  exceeding  0.5  NM  every  35  NM 
can  be  made  at  low  frequencies.  As  a consequence  of  this  effect, 
ray  theory  calculations  of  the  intensity  transmitted  from  point 
to  point  can  show  significant  errors,  even  if  corrections  have 
been  made  for  source  strength. 

An  important  finding  of  this  study  is  that  the  errors  in  ray 
path  have  no  significant  effect  on  the  range-average  of  trans- 
mitted intensity.  Therefore,  they  do  not  affect  the  ambient  noise 
calculated  for  a distribution  of  ships  over  a broad  interval  of 
range.  No  correction  need  for  made  for  this  second  error. 
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1.  TECHNICAL  INTRODUCTION  AND  SUMMARY 

The  accuracy  of  ray  theoretical  calculations  is  founded  in 
large  part  upon  a hypothesis  that  the  fractional  change  of  sound 
speed  in  a distance  equal  to  the  wavelength  is  sufficiently  small. 
Obviously,  the  hypothesis  becomes  increasingly  poor  as  the  sound 
frequency  diminishes.  Furthermore , since  the  largest  gradients 
of  sound  speed  are  observed  near  the  sea  surface  in  typical 
oceanic  environments,  the  errors  will  be  greatest  for  sound  pro- 
pagating near  the  surface.  Here  we  explore  the  nature  of  those 
errors  and  the  possibility  of  correcting  for  them,  giving  parti- 
cular attention  to  the  important  case  where  the  sound  source  is 
near  the  surface. 

The  physical  problem  of  interest  is  that  component  of  low- 
frequency  noise  observed  in  the  deep  ocean  which  is  attributable 
to  surface  ships  distributed  over  the  oceanic  basin.  Surface 
ships  can  create  sound  by  many  different  mechanisms,  some  of  which 
may  be  idealized  as  simple  sources  (monopoles)  and  others  as 
localized  forces  (dipoles)  with  various  orientations.  Some  of 
the  sound  energy  is  refracted  into  the  SOFAR  channel  and  propa- 
gates to  great  distances;  the  rest  interacts  with  the  ocean 
bottom  and  is  attenuated  more  rapidly.  We  shall  ignore  the 
latter  component.  The  noise  observable  at  any  one  point  is 
caused  by  many  ships,  located  at  different  ranges.  Therefore, 
a prediction  of  point-to-point  transmission  is  of  less  interest, 
in  the  end,  than  predictions  valid  for  averages  over  the  range 
variable . 

The  overall  problem  of  long-range  propagation  from  a low- 
frequency,  near-surface  source  includes  two  parts.  The  first 
part  concerns  the  strength  of  radiation  that  is  observed  deep  in 
the  ocean  at  a relatively  short  range.  This  strength  varies 
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with  the  depression  angle  (re  horizontal)  at  which  the  sound  is 
observed  to  be  propagating.  In  ray  theory,  this  strength  can  be 
evaluated  by  adding  coherently  two  components:  one  direct  from 

the  source  and  one  from  the  image  (i.e.,  the  surface-reflected 
path).  A recent  study  by  Pedersen  et  al.  [2_,  indicates  that  this 
strength  may  be  poorly  estimated  by  ray  theory  for  some  realistic 
combinations  of  source  depth,  frequency,  and  sound  speed  gradient. 

The  second  part  of  the  problem  concerns  the  subsequent 
spreading  of  this  energy  as  it  propagates  down  the  channel.  In 
ray  theory,  the  energy  travels  along  ray  paths  whose  location  is 
governed  by  Snell's  law,  and  the  energy  flux  between  neighboring 
rays  is  conserved.  A recent  study  by  Murphy  and  Davis  [2]  indi- 
cates that  the  ray  paths  and  the  distance  between  neighboring 
paths  are  both  distorted  when  low  frequency  sound  passes  through 
regions  having  a large  speed  gradient.  Hence,  ray  theory  may 
lead  to  significant  errors  in  estimating  the  transmission  from 
point  to  point.  However,  we  shall  see  in  the  present  study  that 
this  error  is  largely  compensated  for  when  the  range-averaged 
transmission  is  calculated  by  simple  ray  theory. 

This  study  has  two  objectives.  First,  we  want  to  quantify 
the  domain  of  situations  where  significant  errors  in  shipping 
noise  predictions  can  arise  from  the  use  of  ordinary  ray  theory 
to  calculate  the  sound  generation  and  transmission.  Second,  we 
wish  to  determine  how  corrections,  based  on  wave  theory,  can  be 
introduced  into  ray  theory  calculations  to  compensate  for  those 
errors . 

The  reader  should  particularly  note  two  aspects  of  this 
statement  of  objectives.  First,  it  is  assumed  that  the  test  of 
significance  of  any  errors  is  to  be  found  in  shipping  noise 
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prediction.  The  peculiarities  of  ships  (many  shallow,  low- 
frequency  sources,  widely  distributed  in  range)  strongly  affect 
the  balance  of  attention  given  different  problems.  Second,  it 
is  taken  as  a premise  that  ordinary  ray  theory  calculation  schemes 
will  be  retained,  but  that  corrections  can  and  will  be  made  to  the 
details  of  their  procedures.  Thus,  we  exclude  the  not  entirely 
implausible  view  that,  given  the  errors,  a different  calculation 
scheme  should  be  used  (e.g.,  modal  series). 

At  this  point,  we  summarize  the  analytical  approach  and  the 
principal  results  for  the  technical  reader  who  wishes  to  be  spared 
the  elaboration  of  details  that  follow. 

It  is  assumed  that  the  environmental  parameters  vary  only 
with  depth  and  that  large  speed  gradients  are  limited  to  the 
surface  layers.  We  consider  transmission  from  a shallow  point 
in  the  high-gradient  layers  to  a deep  point  at  long  range. 

The  analytical  formulation  proceeds  by  these  steps,  familiar 
from  Brekhovskikh  [Ref.  4,  Sec.  361  and  many  other  analyses  in 
the  literature: 

i.  The  sound  field  is  expanded  into  plane  waves  inclined 
at  various  angles  to  the  horizontal  at  the  deep 
observation  point. 

ii.  The  separated  wave  equation  governing  the  depth-wise 
distribution  of  sound  pressure  for  a single  plane 
wave  is  "solved",  in  the  sense  that  the  solution  is 
explicitly  defined.  Procedures  for  evaluating  the 
solution  are  put  off  for  later  consideration. 
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ill.  The  responses  to  the  individual  plane  waves  are 

superposed,  yielding  an  integral  expression  for  the 
total  sound  field  at  the  shallow  observation  point. 

At  a later  stage  of  analysis,  reciprocity  will  be 
invoked  to  justify  identifying  the  shallow  observa- 
tion point  as  the  location  of  the  source. 

iv . The  multipath  expansion  of  the  integrand  is  intro- 
duced, so  that  the  total  sound  field  is  expressed 
as  a sum  of  integrals. 

v.  The  integrals  of  the  sum  are  evaluated  by  the  method 
of  stationary  phase  when  conditions  are  such  that  a 
stationary  point  exists;  those  integrals  having  no 
stationary  point  are  dropped. 

Brekhovskikh  demonstrates  that  ordinary  ray  acoustic  result 
from  using  in  step  (ii)  the  ordinary  WKB  expression  for  the 
solution  to  the  separated  Z-dependent  wave  equation.  The  sound 
field  is  found  as  a sum  of  components  identifiable  as  transmissions 
along  the  many  ray  paths  defined  by  Snell's  law. 

The  desired  corrections  to  ordinary  ray  theory  result  from 
using  in  step  (ii)  evaluation  procedures  that  are  more  nearly 
precise  than  ordinary  WKB  procedures.  There  are  various  alterna- 
tives. One  might  use  generalized  WKB  procedures  in  the  manner 
of  Murphy  and  Davis  in  their  development  of  a "modified  ray  theory" 
[Ref.  2].  At  one  level  of  complexity,  this  procedure  amounts  to 
representing  the  sound  speed  profile  as  one  that  is  locally  iso- 
gradient with  slowly  changing  gradient.  The  solution  for  the 
sound  field  is  correspondingly  expressed  in  terms  of  Airy  functions 
(precise  in  a constant  gradient)  with  slowly  changing  parameters. 


A 


4, 


Report  No.  3326 


Bolt  Beranek  and  Newman  Inc. 


Another  alternative  is  to  use  numerical  integration  from  the 
ocean  surface  in  order  to  evaluate  the  solution  to  the  one- 
dimensional, separated  wave  equation. 

In  any  case,  the  requirement  for  a solution  more  nearly 
precise  that  the  ordinary  WKB  expression  can  be  met  by  the  appli- 
cation of  improved  procedures  to  a restricted  zone  of  the  ocean. 
Including  the  high-gradient  surface  layers.  The  gradients  observed 
at  depth  In  the  ocean  are  smaller;  moreover,  the  depression  angle 
of  a ray  originating  at  the  shallow  observation  point  increases 
as  it  reaches  into  the  depths.  Both  these  facts  tend  to  make 
the  ordinary  WKB  expression  more  accurate  at  depth.  The  errors 
of  ordinary  ray  acoustics  in  the  shipping  noise  problem  originate 
when  the  rays  are  near  the  ocean  surface. 

Section  4 below  constitutes  a detailed  examination  of  the 
field  associated  with  a single  multipath  integral  at  the  deep 
observation  point.  It  is  shown  that  the  local  sound  field  at 
depth  follows  the  laws  of  ordinary  ray  acoustics  if  the  ordinary 
WKB  expression  is  locally  accurate.  (It  is  not  assumed  that 
ordinary  WKB  procedures  are  valid  at  all  depths.)  Specifically, 
the  field  lines  of  the  normals  to  phase  fronts  satisfy  Snell's 
law.  The  locus  of  observation  points  having  the  same  stationary 
point  (in  the  stationary-phase  analysis  of  the  integral  over 
radial  wavenumber)  is  a normal  to  phase  fronts  and  satisfies 
Snell's  law.  The  power  in  a ray  bundle  (defined  as  an  interval 
of  radial  wavenumber)  is  invariant  along  such  a locus.  These 
are  all  laws  of  ordinary  ray  acoustics. 

However,  two  errors  are  identified.  First,  the  path  of  a 
ray  (the  locus  defined  above)  is  not  that  predicted  by  ordinary 
ray  acoustics.  Specifically,  an  error  is  made  every  time  that 
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a ray  cycles  through  the  high-gradient  surface  layers,  a result 
carefully  delineated  in  the  prior  work  of  Murphy  and  Davis  [Ref.  2], 
The  range-wise  period  ("cycle  length")  of  such  ray  paths  is  incor- 
rectly evaluated.  Therefore,  the  error  in  path  location  grows 
with  range  in  multicycle  transmission  to  long  ranges.  Moreover, 
since  this  error  in  range  period  varies  from  ray  to  ray,  ordinary 
ray  theory  inaccurately  assesses  the  ray-tube  width,  leading  to 
errors  in  sound  intensity  although  not  affecting  the  power  in  a 
ray  bundle. 

The  second  error  of  ordinary  ray  theory  arises  in  the  cal- 
culation of  sound  power  in  a ray  bundle,  when  and  only  when  an 
observation  point  (e.g.,  the  source)  lies  in  the  high-gradient 
surface  layers.  This  leads  to  more  errors  in  sound  intensity. 

The  second  error  is  very  simply  compensated  for  by  the 
introduction  of  an  effective  source  strength,  defined  as  the 
source  strength  which  yields  the  correct  value  of  sound  power  in 
a ray  bundle  when  used  in  the  formulas  of  ordinary  ray  acoustics. 
These  ideas  are  developed  in  Sec.  5.  The  value  of  the  effective 
source  strength  is  found  from  the  precise  solution  at  shallow 
depths  of  the  separated,  one-dimensional  wave  equation.  Just 
like  Lloyd's  mirror  formula,  it  includes  the  effect  of  the  coherent 
sum  of  what  in  ordinary  ray  acoustics  are  up-going  and  down-going 
rays  with  the  same  depression/elevation  angles  at  the  source. 

The  strength  of  the  sum  is  associated  with  one  down-going  ray. 

The  result  differs  both  from  Lloyd's  mirror  formula  and  from  the 
result  of  strict  ray  tracing  with  coherent  addition. 

Next,  an  examination  is  made  of  the  effect  of  averaging  over 
range  the  transmitted  squared-pressure  -(Sec.  6).  It  will  be 
recalled  that  the  shipping  noise  problem  involves  many  ships. 
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widely  distributed  in  range.  Even  though  calculation  of  point- 
to-point  transmission  may  be  used  as  a numerical  procedure,  the 
desired  answer  involves  an  average  over  many  ranges.  We  have 
adapted  to  the  present  analysis  a range-averaging  procedure 
previously  used  in  the  context  of  ordinary  ray  theory  [Smith, 

Ref.  13].  A statistical  estimate  of  the  contribution  of  each 
ray  bundle  is  formed  on  the  hypothesis  that  every  range  within 
one  range  period  is  equally  probable.  The  resulting  expression 
for  range-av©raged  transmission  is  found  to  be  negligibly  af- 
fected by  the  errors  in  ray  path,  described  above.  However,  it 
is  essential  that  corrections  be  made  for  the  errors  in  the  power 
injected  into  ray  bundles,  by  using  the  "effective  source  strength." 

We  observe  that  the  range-averaged  transmission,  corrected 
by  the  effective  source  strength,  is  easily  calculated  directly 
from  the  analytical  expression  for  it.  It  is  net  necessary  to 
make  many,  erroneous  calculations  of  point-to-point  transmission 
and  then  average  those  results.  Moreover,  this  analytic  averaging 
procedure  could  be  adopted  to  a range-dependent  ocean  by  means 
of  the  adiabatic  invariance  of  action  [cf.  Ref.  13]. 

Numerical  and  analytical  results  for  the  errors  and  cor- 
rection factors  are  developed  in  Sec.  7 for  sources  in  an  iso- 
gradient surface  layer.  The  results  are  based  on  the  precise 
wave  solution  in  terms  of  Airy  functions.  Throughout  the  analysis, 
a depression  angle  computed  from  Snell's  law  is  frequently  used 
as  the  parameter  that  identifies  a ray.  It  is  especially  con- 
venient to  do  so  when  talking  about  corrections  to  ordinary 
ray  theory,  even  though  the  angle  does  not  have  the  same  physical 
meaning  as  in  ray  theory  when  the  wave-theoretical  corrections 
are  significant. 
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The  error  in  the  range  period  of  a ray  varies  with  the 
depression  angle  at  the  surface,  being  a maximum  for  the  ray  that 
is  horizontal  at  the  surface.  Significant  errors  exist  for  rays 
that  vertex  below  the  surface  (within  a wavelength)  as  well  as 
for  those  that  strike  the  surface.  Figures  5 and  6 show  results 
for  typical  values  of  low  frequency  and  high  gradient.  The  maxi- 
mum error,  for  the  grazing  ray,  is  given  by* 

AX  = 1>26 km  , (7^) 

fl/3  g*2/3 

Hz  sec-1 

where  f is  frequency  and  g*  is  the  negative  of  the  sound-speed 
gradient,  hence  a positive  number.  As  discussed  in  detail  in 
Sec.  7>  this  formula  is  not  applicable  for  very  small  values  of 
the  gradient.  However,  within  its  range  of  validity,  errors  in 
excess  of  1 km  are  predicted  for  frequencies  of  200  Hz  or  less 
and  gradients  of  -0.1  sec-1  or  more. 

The  effective  source  strength  of  a unit  monopole  source 
is  found  to  exceed  that  given  by  ordinary  ray  theory  by  6 dB 
and  more  for  small  depression  angles  at  the  source.  For  realistic 
numbers  (f  = 25  Hz,  g*  = 0.5  sec-1,  source  depth  = 20  ft),  the 
error  starts  to  grow  significantly  when  the  depression  angle  is 
reduced  below  8 degrees  (Fig.  8).  Larger  errors  are  found  for 
smaller  values  of  the  ratio  of  frequency  to  gradient  (Fig.  9). 

An  analytical  approximation  for  the  effective  source  strength 
of  a monopole  has  been  found.  For  steeply  inclined  rays,  the 
well-known  Lloyd's  mirror  formula  is  indistinguishable  from 


*The  equation  numbers  are  those  of  the  main  text. 
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ordinary  ray  theory  and  from  more  accurate  results.  Lloyd's 
mirror  formula  for  a shallow  source  whose  strength  is  unity  in 
free  space  yields  a directional  source  strength  given  approximately 
by 

D^M  « 158 ( z/A  ) 2 sin 2 0 , ( 78  ) 

where  z is  source  depth,  A is  the  wavelength,  and  0 is  depression 
angle  at  the  source.  By  contrast,  the  Airy-function  analysis 
yields  an  effective  source  strength  given  by 

D0  « 68(g*/f)l/3  (z/A)2  sin0  . (77) 

This  formula  is  valid  in  the  domain  where  it  exceeds  the  Lloyd's 
mirror  result,  that  is  for 

g*/f  > 12.5  sin 3 0 . (79) 

This  is  the  domain  where  both  ordinary  ray  theory  and  Lloyd's 
mirror  formula  are  significantly  in  error. 

Similar  errors  of  ordinary  ray  theory  are  found  in  the  case 
of  dipole  (force)  sources.  The  details,  which  vary  according  to 
the  orientation  of  the  force  (vertical,  horizontal,  etc.),  are 
developed  in  Sec.  7.  The  domain  in  which  wave-theoretical  cor- 
rections to  the  effective  source  strength  must  be  used  is 
approximately  the  same  as  Eq . 79  for  the  monopole  source.  More- 
over, the  correction  factor  that  must  be  applied  to  the  Lloyd's 
mirror  formula  for  a dipole  is  approximately  the  same  as  in  the 
monopole  problem. 
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2.  INTEGRAL  FORMULATION 

We  are  concerned  with  sound  propagation  from  pure-tone 
volume  and  force  sources  (monopoles  and  dipoles)  located  in  an 
inhomogeneous  fluid  whose  sound  speed  c(z)  varies  with  the  depth 
coordinate  z and  whose  density  p is  constant  [3].  All  properties 
are  independent  of  horizontal  position.  Using  complex  variables, 
we  include  the  time  dependence  implicitly  as  a factor  exp(-iwt), 
where  m = 2iTf  is  the  frequency. 

The  sound  pressure  p(R)  and  particle  velocity  vector  u(R), 
observable  at  a vector  position  R in  response  to  a volume  source 
Q and  a vector  force  source  F both  located  at  R = R0,  must  satisfy 
the  equations  of  continuity  and  of  force  balance: 

-iujc_2p  + pV*u  = pQ6(R-RQ)  (1) 

-iwpu  + Vp  = F6(R-Rq)  , (2) 

where  6(***)  indicates  a three  dimensional  delta  function.  The 
particle  velocity  is  readily  eliminated,  with  an  inhomogeneous 
wave  equation  for  the  pressure  as  the  result: 

V2p  + k2p  = iojpQ6(R-R0)  - F • V Q 6 ( R-R  0 ) , (3) 

where  k(z)  = oj/c(z).  We  have  used  the  vector  identities 
V • [F6 ( R-R  Q ) ] = F*V6(R-Ro)  = -F*V05(R-R0)  , 

where  F is  a constant  vector  and  V0  denotes  the  gradient  with 
respect  to  the  source  location  Rq.  The  complete  mathematical 
statement  of  the  problem  requires  additional  specification  of  the 
boundary  conditions;  see  below. 
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The  field  in  response  to  either  source  is  readily  derived 
from  the  Green's  function,  defined  as  the  solution  to 

[V2  + k2  ( z ) ] G ( R , R ) = -4tt6(  R-R0)  , (4) 

with  the  same  boundary  conditions.  This  Green's  function  ap- 
proaches | R— R 0 | 1 as  | R-R0  | -+  0,  and  thus  has  the  standard  normal- 

ization to  the  field  at  a unit  range.  Under  broad  conditions, 
which  we  assume  to  hold  here,  G is  reciprocal,  in  the  sense  that 
its  value  is  invariant  to  an  interchange  of  source  and  receiver 
locations . 

Cylindrical  coordinates  (r,z,a)  are  convenient  for  further 
analysis.  We  take  the  axis  r=0  to  be  a vertical  line  through  the 
source.  The  origin  of  depth  is  placed  at  the  sea  surface,  with  z 
increasing  downward.  The  source  depth  is  z1.  The  solution  for  G 
is  independent  of  a.  A typical  sound  speed  profile  is  sketched 
in  Fig.  1(a).  In  this  analysis,  the  source  is  deep;  we  shall 
seek  solutions  for  points  near  the  surface,  relying  on  reciprocity 
to  justify  using  the  solution  for  the  actual  problem  of  interest 
with  the  source  near  the  surface.  The  sea  surface  is  taken  as  a 
boundary  at  which  G=0. 

The  boundary  condition  at  the  bottom  will  be  left  vague.  We 
may  anticipate  that  it  plays  only  a minor  role  in  long-range 
propagation  since  we  are  interested  only  in  ducted  energy:  i.e., 
rays  that  do  not  touch  the  bottom,  or  trapped  modes.  Bottom 
properties  must  be  specified  in  order  numerically  to  perform  a 
modal  analysis,  but  to  the  small  extent  that  they  affect  results 
the  effect  is  probably  erroneous  due  to  artificiality  of  the 
modeling  of  the  bottom. 
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(A):  ACTUAL  SOUND  SPEED  PROFILE.  (B)  AND  (C):  SEMI-INFINITE 

PROFILES  USED  IN  AUXILIARY  PROBLEMS  TO  DEFINE  REFLECTION 
COEFFICIENTS  FOR  THE  UPPER  AND  LOWER  SPACES,  RESPECTIVELY 
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Brekhovskikh  has  shown  that  the  solution  for  G can  be 

expressed  as  a superposition  of  cylindrical  waves  diverging  from 

the  source  [6].  Each  wave  is  characterized  by  its  value  of  the 

radial  wavenumber  k , and  the  formal  solution  is  an  integral  over 

all  real  values  of  k : 

r 


G(z2,r,z1 ) = ^ 


l’u,t(z2>zi>kr) 


1 + V„  (z  ,k  ) 

i ,u  i ’ r 

1 - V„V 
£ u 


( i ) % 

IT  ' (k  r)  7 — i- — v dk  , 
o r k ( z ) r ’ 


(5) 


where 


kz ( z)  = [k2 (z  ) - k£] 2 . (6) 

The  first  subscript  of  the  pairs  appearing  in  the  integrand  is 
used  when  z2 < z1  and  the  second  for  z2  > . 

The  terms  appearing  in  brackets  in  Eq.  5 (g  are 

found  from  solution  of  two  auxiliary  problems:  the  reflection 

of  an  incident  plane  wave  of  unit  amplitude  from  either  the 

upper  (u)  or  lower  (£)  half  of  the  channel  ( z £ z j or  z >_  z j ) . 

Consider  the  reflection  analysis  of  the  upper  half.  The  actual 

channel  for  z < zx  is  supplemented  by  a semi-infinite  homogeneous 

halfspace,  Zj  < z < with  constant  sound  speed  c(z1).  The 

resulting  sound  speed  profile  is  shown  in  Fig.  1(b).  A plane 

wave  is  generated  in  the  homogeneous  half space,  traveling  upward 

at  an  angle  0 with  respect  to  horizontal,  where  tan9  = k (z  )/k  . 

z i r 

This  incident  wave  has  unit  amplitude  and  is  phased  to  have  zero 
phase  angle  at  z,.  Reflection  of  the  incident  wave  will  cause 
response  in  the  channel  and  a reflected  wave  traveling  at  an 
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angle  -0  In  the  homogeneous  half space.  Thus,  the  total  pressure 
field  has  the  form 

-ik  (z-z  ) ik  (z-z  ) 

h 1 zi  1 

e + V e 

u 

ik  x 

Su(z)  e r , z < zi 
where  we  write  k = k (z.).  Hence,  V is  the  complex  plane- 

Z j Z 1 \X 

wave  reflection  coefficient  for  the  upper  half  channel  and 
is  the  corresponding  response  inside  the  half  channel. 

The  terms  g^  and  are  defined  by  a similar  reflection 
problem  in  which  the  lower  half  of  the  actual  channel  is  supple- 
mented by  a semi-infinite  homogeneous  half space  above  it.  This 
sound  speed  profile  is  shown  in  Fig.  1(c). 

A mathematician  might  take  another,  totally  equivalent  view. 
In  a straightforward  derivation  of  Eq . 5,  the  Hankel  transform  of 
Eq.  4 is  taken  to  eliminate  the  r-dependence , yielding  an  ordinary 
differential  equation  in  z with  a source  term  -25(z-Zj).  Its 
solution,  which  is  the  bracketed  term  in  Eq.  5>  can  be  constructed 
from  the  solutions  to  boundary-value  (initial-value)  problems  for 
the  upper  and  lower  halves  of  the  channel.  For  the  upper  half  we 
require  the  solution  to 

g"  + k2  g = 0 , k2  = k2(z)  - k2  , (8) 

&u  z &u  z r 

with  appropriate  boundary  conditions  at  z=0.  (The  prime  denotes 
differentiation  with  respect  to  z.)  The  initial  value  at  z=z^ 
is  defined  as 


ik  x 
r 


z > z 


,1 

I 


(7) 
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1 + V 

u 


V 

u 


t(-ikz  - 


«yKu)/(-ikz 


Ei/gu)]z-z 


(9) 


The  functions  g^,  are  found  from  a similar  initial-value  prob- 
lem for  the  lower  half  of  the  channel,  except  that  in  the  defini- 
tion of  the  reflection  coefficient,  Eq . 9,  ikz  is  substituted  for 
-ik^.  The  desired  solution  to  the  inhomogeneous  equation  (with 
source  term)  is  represented  by  g and  g^  when  they  are  scaled  to 
satisfy  the  appropriate  continuity  conditions  across  z=z 
(continuous  function  and  discontinuous  first  derivative).  Then 
the  inverse  Hankel  transform  leads  to  Eq . 5>  when  symmetry  of  the 
integrand  is  invoked  to  extend  the  lower  limit  of  the  integral 
to  -°°  (cf.  Ref.  4,  p.  243). 
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3.  MULTIPATH  EXPANSION 

Brekhovskikh  has  shown  how  the  general  Integral  expression 
for  the  sound  field,  Eq . 5,  is  readily  expanded  into  a sum  of 
terms,  each  of  which  can  be  shown,  with  appropriate  approximations, 
to  represent  the  contribution  of  energy  travelling  along  a distinct 
ray  path  from  the  source  [£].  In  the  present  development  we 
follow  Brekhovskikh,  with  minor  generalizations. 

The  essential  step  toward  a multipath  expansion  is  to  expand 
the  factor  in  the  brackets  of  Eq.  5 into  a power  series  in  the 
reflection  coefficients,  using  the  identity 


[1  - v v r1  - [ <vuvn  . no) 

n=0 

which  is  convergent  in  any  physical  case  since  slight  dissipation 
causes  |v|  to  be  absolutely  less  than  unity.  The  Green's  func- 
tion then  has  the  form 

oo 

G = l [G+  + G~]  , C 11 ) 

n=0  n n 


where  (for  z2  < Zj)  the  typical  term  is 


Gn(z2,r,z,)  = (i/2) 


gu  v^vj  H['>(krr)  [krAz(z,>]dkr 


(12) 


The  integrand  for  G differs  in  having  VQ  raised  to  the  power 

n ~ 

n+l.  The  expressions  for  z2  > z differ  in  having  the  subscripts 
u and  H interchanged. 
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r 


Rrekhovskikh  shows,  for  a special  case,  that  the  term  G+ 

n 

corresponds  to  transmission  along  the  ray  path  that  leaves  the 
source  in  an  upward  direction  and  loops  n times  into  both  the 
upper  and  lower  halves  of  the  channel,  being  turned  each  time  by 
reflection  or  refraction.  The  path  corresponding  to  G~  leaves 
the  source  in  a downward  direction,  being  turned  n+1  times  in  the 
lower  half  and  n times  in  the  upper  half  channel.  (See  Fig.  2.) 

In  the  case  Rrekhovskikh  chooses  for  illustration  of  the 
technique,  the  boundary  z=0  is  rigid  and  the  receiver  lies  on 
the  boundary.  However,  the  generalization  is  apparent.  Within 
the  WKR  approximation,  the  field  at  the  receiver  point  can  be 
split  into  upward  and  downward  travelling  parts,  whose  inter- 
ference is  described  by  the  function  g . If  G+  were  correspond- 
ingly split  into  two  integrals,  each  would  be  found  to  describe 
a different  ray  path.  One  of  the  paths  would  approach  the  receiver 
from  below  and  the  other  from  above,  after  reflection  from  the 
surface  (see  insert  in  Fig.  2). 

However,  in  the  cases  of  present  interest,  where  the  receiver 
is  near  the  surface  and  at  a long  range,  it  would  be  inappropriate 
to  attempt  such  a decomposition  of  G . The  two  parts  are  not 
resolved  in  any  measurements  of  interest.  The  analysis  must 
retain  the  coherent  interaction  between  the  parts,  as  in  Eq.  12. 

The  correspondence  between  the  integral  for  Gn  and  the  pres- 
sure transmitted  along  the  ray  path  is  established  by  (1)  using 
expressions  for  g , V , V.  based  on  the  WKB  approximation  to  solu- 

U U X/ 

tions  of  Eq . 8,  and  (2)  evaluating  the  integral  by  the  method  of 
steepest  descents  in  the  complex  k^  plane.  The  resulting  values 
of  amplitude  and  phase  are  found  to  be  the  same  as  in  ordinary 
ray  theory,  except  foi*  a negligible  correction  to  the  phase. 
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In  the  present  study  we  hypothesize  the  failure  of  the 
of  the  ordinary  WKB  approximation  in  the  near-surface  layers 
(where  the  receiver  point  is  located),  as  a consequence  of  a low 
frequency  and  large  speed  gradient.  However,  we  hypothesize  the 
validity  of  the  WKB  approximation  near  the  source  (z=z, ) and, 
indeed,  throughout  the  depth  of  the  ocean  except  the  surface 
layers.  We  follow  Brekhovskikh’ s procedure,  insofar  as  possible, 
and  develop  corrections  to  ordinary  ray  theory. 

First  let  us  develop  the  formal  expression  for  G , Eq . 12, 
using  the  method  of  stationary  phase.  We  can  justify  a posteriori 
using  the  asymptotic  form  for  the  Hankel  function: 

hJ°(z)  ~ (2/ttz)1'5  ei(z_1T/i0  . (13) 

The  phases  of  other  terms  in  the  integrand  are  denoted  by: 
arg  VUV£  = <p(  zi  ,kr) 

arg  gu  = y(z2 ,Zj ,kr)  . (14) 


[In  the  case  of  a term  G~,  y would  denote  arg(g  V^).]  We  take  as 
the  value  of  the  integral  that  given  by  the  method  of  stationary 
phase: 


A(x)e1<J>^x^dx 


A( x0  )e 


i [$  ( x 0 )±tt/2] 


where  xQ  is  the  point  of  stationary  phase  where  the  derivative 
$'  vanishes,  and  the  choice  of  sign  in  the  exponent  agrees  with 
the  sign  of  <i>"  [2],  The  total  phase  of  the  integrand  in  Eq . 12 
is 


$ = k r + n<J>  + y - tt/4 


(16) 
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which  has  a stationary  point  at  the  root,  denoted  by  ky,  of  the 
equation 


0 = r + n<t>’  + y'  at  kr  = kv  , (17) 

where  the  prime  indicates  partial  differentiation  with  respect 

to  k . 
r 


The  straightforward  combination  of  these  relations  yields 
for  the  squared  magnitude  of  the  multipath  term: 


Gn(z2’r’zi) 


2 _ 


g ( z , z )|2  | V V 0 j n k 
Bu  2 1 1 1 U X.  1 I 


r k2  ( z ) I $>" 
z i ' 


(18) 


k =k 
r v 


and  for  its  phase: 


arg  G = <t>,  _R  + it  + tt/2  , 
r v 


(19) 


with  the  sign  chosen  to  agree  with  the  sign  of 


i 

I 

\ 

] 
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4.  LOCAL  RAY  BEHAVIOR 

In  this  section  we  examine  the  spatial  dependence  of  the 
proceeding  expressions  for  the  contribution  of  one  multipath 
term  in  order  to  find  the  similarities  and  differences  with  the 
results  of  ordinary  ray  acoustics.  The  physical  discussion  is 
simpler  if  we  once  more  appeal  to  reciprocity:  we  shall  call  the 

near-surface  point  at  z2  the  "source",  and  the  deep  point  at  z 
the  "receiver". 

The  results  of  this  section  are  intimately  related  to  the 
validity  of  the  WKB  approximation  to  the  solutions  g^  and  g^  of 
the  z-dependent , one-dimensional  wave  equation  (Eq.  8).  The  local 
variations  of  sound  field  may  be  well  described  by  the  WKB  solu- 
tion in  some  intervals  of  z,  which  we  call  "WKB  regions".  The 
WKB  solution  may  be  inaccurate  in  other  intervals,  called  "non- 
WKB  regions".  In  what  follows,  we  assume  that  the  surface  layers, 
containing  the  source  at  z2,  constitute  a non-WKB  region,  while 
the  rest  of  the  ocean,  containing  the  receiver  at  z , is  a WKB 
region.  The  results  could  readily  be  generalized  to  cases  with 
several  non-WKB  regions  [fl]. 

The  WKB  approximation  represents  the  general  solution  to  the 
one-dimensional  wave  equation  as  the  sum  of  components  travelling 
upward  and  downward,  in  the  form 


g(z) 


[kz(z  )]"** 


z 

k dz 
z 

+ 


+i 

C e 
2 


(20) 


where  Cj  and  C2  are  arbitrary  constants  [9].  In  the  case  of  the 

function  gu  (Eq.  8),  the  constant  Cj  is  by  definition  to  be 

chosen  so  that  the  upgoing  wave  is  unity  at  z=z,.  The  constant 

C is  related  to  the  reflection  coefficient.  Then  the  solution 
2 

can  be  written 
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Ru(.z,z1  ) 
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kz(z,) 
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k dz 
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+ Vu(zi)e 


k dz 
z 


(21) 

for  any  two  points  z and  z,  both  in  the  WKB  region.  If  the  WKB 

region  were  to  extend  all  the  way  to  the  surface  z=0  and  k2(0)>0, 

z 

then  it  can  easily  be  shown  that 


rzi 

= -exp 

i2 

k dz 

0 Z 

g ( z , z , ) = -i2sin 
u i 


k dz 
z 


kz(z!)/kz(z) 


l-5 

rzi 

exp 

i 

k dz 
z 

0 

(22a) 

This  corresponds  to  a ray  that  reflects  from  the  surface.  On  the 
other  hand,  if  k2(0)<0  corresponding  to  a ray  that  vertexes  below 
the  surface,  then  the  ordinary  WKB  expressions  are 

V (z,  ) = 
u i 


gu(z,z,) 


= -expi2 

71  + 

If  + J 

fz»  , 

kzdz 

2 

7 

Jv 

= i2  sin 

77  + 
2f  + 

kzdz| 

kz(zi )/kz(z) 

2 

J 1 

expi 


71 

r + 

kzdz 

' z ' 

(22b) 


where  zy  is  the  depth  of  the  turning  point,  defined  by  kz(zv)=0  [20] 


22 


1 

Bolt  Beranek  and  Newman  Inc. 


Now  we  wish  to  find  the  dependence  upon  zx  of  g (z2,Zj) 
when  z.  is  in  the  WKB  region  but  z,  is  not.  The  function  g is 
by  definition  the  response  of  a fixed  system  to  a unit  upgoing 
wave  at  zx;  its  dependence  upon  z2  arises  solely  from  the 
necessity  to  renormalize  to  a unit  upgoing  wave  when  zx  is 
changed.  Therefore  the  dependence  of  g upon  z,  is  completely 
given  by  the  factor  in  braces  in  Eq . 21,  which  describes  the 
upgoing  wave  in  the  WKB  region.  In  particular,  the  phase 
Y=arg  gu  satisfies 

3y(z,zi ,kr)/3zi  = kz(zi)  , (23) 

and  the  magnitude  satisfies 

3[ |gul 2/kz(z1 ) ]/3Zj  = 0 . (24) 

But,  by  the  same  argument,  the  last  factor  of  Eq . 21  must  be  inde- 
pendent of  Zj . Therefore,  the  reflection  coefficient  satisfies 

9|Vj/9zi  = 0 , 3arg  Vu/9Zi  = 2 k^Zj)  . (25) 

When  the  same  argument  is  carried  out  for  the  lower  half  channel, 
the  signs  in  both  exponents  of  Eq . 21  must  be  changed.  It  follows 
that 

9|Vjl|/9z|  = 0 , 9arg  V£/3z]  = -2  kz ( z , ) , 

9 1 V V.|/9z  = 0 , 9<J>/9z  = 0 where  4)  = arg  V V. 

(26) 

We  are  now  in  a position  to  examine  the  variations  with 
receiver  location  (coordinates  r and  z ) of  the  expressions  given 
in  Eqs.  18  and  19  for  the  strength  of  a single  multipath  term. 


A 


Report  No.  3326 


w 


Report  No.  3326 


Bolt  Beranek  and  Newman  Inc. 


Many  similarities  and  some  differences  are  found  in  comparison 
with  the  corresponding  results  of  ordinary  ray  acoustics. 

Consider  the  locus  of  positions  having  the  same  value  k 

v 

of  the  radial  wavenumber  at  the  stationary  phase  point.  Equation 

17  determines  k as  a function  of  r and  other  variables:  its  in- 
v 

version  yields  an  expression  for  r as  a function  of  kv : 

r(Zj >z2  >kv’n)  = -3[n4>(l<v)  + y(Zj  ,z2  ,kv)]/3kv  . (27) 

The  partial  derivative  with  respect  to  is  readily  found,  using 
Eq.  23: 


Br/SZj  = -92Y/3z13kv  = -3kz ( z l , ky) /9kv  = kv/kz(Zj)  , 

(28) 

since  k2  = k2  - k2  by  definition.  But  this  is  just  the  equation 
for  an  ordinary  ray  path  governed  by  Snell's  law  with  a vertexing 
sound  speed 


cv  = m/kv  . (29) 

The  tangent  to  the  locus  is  depressed  from  the  horizontal  by  an 
angle  0j(z  ) that  satisfies 

kv  = k(z l ) cose  j , 

k = k( z j ) sin0  j . ( 30 ) 

Moreover  the  field  lines  of  the  gradient  of  phase  (Eq.  19) 
describe  the  same  locus,  for  it  is  readily  shown  from  Eqs.  16, 
23,  and  26  that 

9<1 >/9r  = kv  , 3$/9z1  = ay/SZj  = kz(z1)  . (31) 
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However,  this  ray-like  behavior  is  valid  only  for  local 
variations  of  z i within  the  WKB  region.  Equations  23  and  26  for 
the  z -dependence  of  the  phases  y and  4>  both  fail  in  a non-WKB 
region.  Equation  27  indicates  that  the  locus  of  constant  kv  is 
periodic  in  range  with  a period 

Ar/An  = -3<fr(kv)/3kv  & X(ky)  , (32) 

but  this  range  period  ("cycle  length"  or  "bounce  distance")  is 
not  correctly  predicted  by  WKB  analysis.  Moreover,  the  range 
found  by  ordinary  ray  (WKB)  analysis  has  an  additional  error  due 
to  the  non-WKB  behavior  of  y.  The  locus  of  constant  k^  and  the 
ordinary  ray  path  with  c chosen  as  in  Eq . 29  are  sketched  in 
Fig.  3-  The  error  increases  with  increasing  range  because  of  the 
error  in  the  range  period  X.  The  sketch  exaggerates  the  effect. 
Calculations  for  typical  cases  made  by  Murphy  and  Davis  [72] 
show  errors  of  the  order  of  1 km  or  less  in  a range  period  X 
which,  in  the  deep  ocean,  is  typically  about  70  km. 

In  summary,  we  have  found  that,  within  the  WKB  region,  the 
locus  of  constant  kv  is  the  normal  to  the  phase  fronts,  and  is 
the  same  as  the  path  of  the  ordinary  ray  having  vertex  speed 
cv  = w/kv  except  for  a range-wise  displacement.  The  error  in- 
creases with  range  because  the  range  period  X differs  between 
the  two  paths. 

Consider  next  the  spatial  variations  within  the  WKB  region 
of  the  squared  magnitude  of  pressure,  as  given  by  Eq.  18.  We 
first  calculate  the  power  associated  with  the  sound  field  in  a 
ray  bundle,  where  "ray"  is  interpreted  as  a locus  of  constant  k . 
Because  these  loci  are  equivalent  to  ordinary  ray  paths  in  the 
WKB  region,  the  calculation  procedure  parallels  that  of  ordinary 
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SOURCE 


FIG.  3.  SKETCH  OF  PATHS  OF  MODIFIED  RAY  (LOCUS  OF  CONSTANT  kv)  AND 

CORRESPONDING  ORDINARY  RAY.  THE  MODIFIED  RAY  PATH  IS  NOT  DEFINED 
IN  THE  NON-WKB  REGION  NEAR  THE  SURFACE. 
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ray  theory,  where  the  power  in  a ray  bundle  is  found  to  be  con- 
stant. The  power  P satisfies 

9P/9kv  = I(kv)  9A/9ky  , (33) 

where  I is  the  acoustic  intensity  and  A is  the  transverse  area 
of  a ray  bundle  6 kv  measured  in  a direction  normal  to  the  rays. 
The  intensity  is  related  to  the  squared  pressure: 

I = |Gn|2/pc(z1)  = k(z1)|Gn|2/pm  . (3*0 

(This  formula  yields  the  peak  intensity,  not  the  time-averaged 
value  which  is  half  as  large.  Throughout  this  report  we  use 
peak  intensity,  peak  power,  and  peak  squared  pressure.  The 
formulas  relating  them  are,  of  course,  identical  to  the  formulas 
relating  the  time-averaged  quantities.)  The  transverse  area  is 
found  by  projection  of  the  range  derivative  onto  the  normal  to 
the  rays: 


9A/9kv  = 2iTr  | 9r/9kv  | sin6  ^ 

= 2Trr[kz(z1  )/k(Zj  ) ] |9r/9kv|  , (35) 

where  the  factor  2ffr  arises  from  the  cylindrical  symmetry  of  the 
sound  field  about  a vertical  axis  through  the  source. 

Consistently  with  our  interest  in  ducted  ray  paths,  we 
assume 


iv*  I - 1 


(36) 
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It  follows  directly  from  Eqs . 16  and  27  that  the  range  derivative 
is 


3r/3kv|  = 4>"(kv) 


(37) 


which  is  a factor  in  | G | 2 . The  straightforward  combination  of 
these  results  yields 


3P 

3k 

v 


2tt 

pw 


I S,,  C : 


i»kv}l 


k ( z ) 
z i 


(38) 


where  z2  is  source  depth  and  z is  receiver  depth.  However,  we 
have  already  found  (Eq.  24)  that  this  expression  is  independent 
of  Zj  so  long  as  Zj  is  in  the  WKB  region.  In  brief,  the  energy 
flux  in  a bundle  of  modified  rays  is  invariant  with  position  [12], 
just  as  in  ordinary  ray  acoustics. 


However,  ordinary  ray  acoustics  does  not  correctly  predict 
this  power.  The  factor  I Su I 2 is  not  correctly  predicted  when  the 
source  at  depth  z2  lies  in  the  non-WKB  region.  The  error  is 
independent  of  receiver  depth  z,  (see  Eq . 24)  and  independent  of 
range  (it  is  the  same  for  each  multipath  term,  independent  of  n). 


Ordinary  ray  theory  will  lead  to  other  serious  errors  in  the 
prediction  of  squared  pressure  transmitted  to  a specific  receiver 
point.  First,  the  wrong  ray  (the  wrong  value  of  kv)  will  be 
identified  as  connecting  two  points,  because  of  the  errors  in  ray 
path  elucidated  above.  Second,  since  the  errors  in  ray  path  vary 
as  k changes  from  path  to  path,  the  wrong  value  of  the  range 
derivative  will  be  calculated  from  ordinary  ray  acoustics.  Both 
errors  affect  the  magnitude  of  pressure. 
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5.  EFFECTIVE  DIRECTIONAL  SOURCE  STRENGTH 

In  the  last  section,  we  showed  that  the  sound  power  in  a ray 
bundle  is  invariant  with  depth  and  range  if  the  observation  point 
is  in  the  WKB  region.  (Dissipative  processes  are  ignored  in  this 
analysis.)  However,  we  found  that  ordinary  ray  theory  does  not 
correctly  evaluate  that  power  when  the  source  is  in  a non-WKB 
region  like  the  high-gradient  surface  layer. 

In  this  section,  we  develop  expressions  for  an  effective 
(directional ) source  strength , defined  as  the  source  strength 
which  yields  the  correct  value  of  sound  power  in  a ray  bundle 
when  the  spreading  is  calculated  by  ordinary  ray  acoustics.  Thus, 
the  modification  corrects  for  the  errors  in  ordinary  ray  acoustics. 
The  effective  source  strength  includes  the  coherent  sum  of  what 
in  ordinary  ray  theory  would  be  the  direct  and  the  surface  re- 
flected (or  refracted)  rays.  This  total  strength  is  associated 
with  a single,  down-going  ray  at  the  source. 

Different  expressions  are  required  for  monopole  sources 
(fluctuating  volumes)  and  for  dipole  sources  (fluctuating  forces). 

Ray  Bundle  Power  in  Ordinary  Ray  Acoustics 

Consider  a monopole  source  with  source  strength  S,  where 
"source  strength"  has  the  conventional  meaning  of  the  squared 
pressure  observed  at  a unit  range  from  the  source  in  an  infinite 
medium.  Vie  wish  to  find  the  expression  for  power  in  a ray  bundle 
according  to  ordinary  ray  acoustics.  Since  this  power  is  every- 
where invariant  with  range,  the  calculation  is  most  easily  per- 
formed at  a unit  range  near  the  source,  where  the  sound  is 
spreading  spherically. 
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For  the  source  of  strength  S located  at  depth  z2 , the  in- 
tensity at  a unit  range  is  just 

I = S/pc ( z 2 ) = S k(z2)/po)  , (39) 

independent  of  the  ray  direction.  Now  consider  a ray  bundle 
where  the  rays  are  identified  by  the  ray  invariant 

kv  = w/cv 

where  c v is  the  vertexing  sound  speed.  The  relationship  to  the 
local  depression  angle  6 of  the  ray  is  given  by  Snell's  law 

kv  = k(z)  cos0(z)  ; (40a) 

we  shall  also  use  the  vertical  component  of  wavenumber  defined  by 
kz(z)  = k(z)  sin0(z)  . (40b) 

The  differential  of  power  dP  in  a ray  bundle  dkv  satisfies 
dP/dkv  = I dA/dkv 

where  dA  is  the  corresponding  differential  of  area  on  the  unit 
sphere : 


dA  = 2 tt  cos0  d0  = ( 2 ttt  cos0/k  sin0)dkv 

The  combination  of  these  results  with  Eq . 39  yields  the  desired 
result : 


dP  = 2tt  kvS 
dkv  p<D  kz(z2 ) 


(41) 
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This  expression  gives  the  power  in  a ray  bundle  radiated  from  a 
source  of  strength  S,  according  to  ordinary  ray  acoustics. 

We  now  compare  the  results  of  the  earlier  analysis  which 
includes  the  effects  of  both  surface  interaction  ("reflection") 
and  wave-theoretical  corrections. 

Monopole  Source 

Consider  a monopole  source  whose  intrinsic  source  strength  is 
SQ.  By  this  it  is  meant  that  a squared  pressure  equal  to  S0  would 
be  observed  at  a unit  range  if  the  source  were  located  in  an  in- 
finite homogeneous  medium.  The  Green's  function  G of  Eq.  4 is 
normalized  so  as  to  have  S0  = 1.  On  the  other  hand,  a monopole 
whose  volume  velocity  is  Q has  an  intrinsic  source  strength 


S0  = | mpQ/4ir  | 2 , (42) 

as  is  readily  seen  by  comparison  of  Eqs.  1 and  3. 

The  power  radiated  into  a ray  bundle  by  a monopole  of  unit 
source  strength  was  derived  in  the  preceding  section,  Eqs.  33 
through  38.  For  a monopole  of  strength  S , the  result  is 


9F  ^ 2ttS  o ky(gu(z2  ,Zj  ,kv)  | 2 
9k^  ~ “pm  kz(Zj ) 


(43) 


where  z2  is  source  depth  and  the  expression  has  been  shown  to  be 
independent  of  receiver  depth,  z1.  A comparison  of  Eq . 43  with 
the  corresponding  expression  for  ordinary  ray  acoustics,  Eq.  4l, 
shows  that  the  effective  source  strength  of  the  monopole  is 
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S 


eff 


S D 


o o 


y 


D = 


|SU(Z2 


>kv)  I 


kz(z2)/kz(zx) 


(44) 


When  the  environment  is  such  that  the  WKB  solution  is  accurate 
throughout  the  ocean,  the  directional  factor  DQ  reduces  to  familiar 
expressions  of  ordinary  ray  theory.  From  the  WKB  expressions  for 
gu  given  in  Eq . 22,  we  find 


D 


WKB 

o 


4 sin2 


kzdz 


when  kz2(0)>0,  and 


D 


WKB 

o 


fZ, 

= 4 sin2 

11  + 
Tf  + 

kzdz 

(45a) 


(45b) 


when  k2(zv)  = 0 for  some  z^  > 0 and  < z2.  These  expressions  give 
the  ray-acoustical  modification  to  source  strength  consequent  to 
taking  the  coherent  sum  of  the  contributions  of  the  two  rays 
that  leave  the  source  at  angles  +0.  In  Eq.  45a,  the  upgoing  ray 
reflects  from  the  surface;  in  Eq . 45b,  the  upgoing  ray  is  re- 
fracted without  reaching  the  surface. 


These  ray-acoustical  results  are  to  be  compared  with  the 
familiar  Lloyd's  mirror  formula , which  is  the  precise  solution 
in  an  environment  having  constant  sound  speed: 

Df  = 4 sin2  [kz(z2  )z2  ] . (46) 
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Horizontal  Dipole 

Consider  now  the  radiation  from  a horizontal  dipole  source  - 
that  is,  a fluctuating  force  F,  oriented  horizontally. 

It  is  immediately  evident  from  a comparison  of  the  wave 
equations,  Eqs.  3 and  4,  that  the  field  radiated  from  the  dipole 
is  related  to  the  unit  monopole  field  G by 

„ _ F cosa  9G  ,hr7\ 

P W-  9F  ’ (^7) 

where  a Is  the  azimuth  angle  in  the  horizontal  plane  between  the 
force  vector  and  the  line  from  the  source  to  the  receiver.  But, 
the  dominant  dependence  of  G on  r,  in  the  far  field,  is  in  the 
phase  (Eqs.  16  and  19): 

G “ exp  ikvr 


Hence,  we  have 


ikyF  cosa 


Now,  the  analysis  for  the  power  in  a ray  bundle  exactly 
parallels  the  earlier  analysis  (Eqs.  33  to  38)  for  a unit  monopole 
source.  The  result  for  the  horizontal  dipole  is  readily  seen  to 
equal  the  result  for  the  unit  monopole  multiplied  by  the  squared 
magnitude  of  the  factor  in  brackets  in  Eq . 48.  The  effective 
source  strengths  are  similarly  related.  Thus,  the  effective 
source  strength  of  the  horizontal  force  can  be  written  as 

Seff  = Ck(z2  ) | F| /4tt]2  Dpp  , 


Dh„  = cos20(z  ) cos2a  D0  , 
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DQ  being  given  by  Eq.  144.  The  ray-acoustical  approximation  for 
D0  , Eq . 45,  and  Lloyd's  mirror  formula,  Eq . 46,  are  equally 
applicable  here. 

Vertical  Dipole 

Consider  now  a vertical  dipole  (force)  source  of  strength  F. 
The  comparison  of  Eqs.  3 and  4 shows  that  the  radiated  pressure 
is  related  to  that  for  a unit  monopole  by 

p = [F/4it]  3G/dz2  . (50) 

The  analysis  for  power  in  a ray  bundle  parallels  that  for  a unit 
monopole.  The  result  is  readily  seen  to  differ  only  by  the 
squared  magnitude  of  the  factor  F/4u  and  by  the  substitution  of 
(3gu/3z2)  for  . Thus,  the  effective  source  strength  of  the 
vertical  force  can  be  written  as 


eff 


D. 


VF 


[k(z  2 ) |F’|/4tt]2  Dvp 

kz(z2)  I 3 Su  C z 2 , z ! , k v ) / 3 ssj 
k (z  j ) k2 (z  ) 


(51) 


The  ray-acoustical  approximation  for  Dyp  is  found  by  dif- 
ferentiating the  WKB  expression  for  g , Eq . 22.  For  the  case 
of  rays  that  strike  the  surface  [k2(0)>0],  the  result  is  found 
to  be 


D 


WKB 

VF 


4 sin2 0 ( z2  ) cos 


(52) 
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where  we  have  omitted  a term  involving  (9kz/3z2)  that  is  only 
important  when  the  cosine  is  nearly  zero.  In  the  case  analogous 
to  Lloyd's  mirror,  when  kz  is  independent  of  depth  in  the  layer, 
the  ray-acoustical  result  Is 

Dyp  = 4 sin20(z2)  cos2 [kz ( z2 ) z2  ] . (53) 


r 
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6.  RANGE-AVERAGED  MODIFIED  RAY  THEORY 

Some  of  the  errors  of  ordinary  ray  theory  have  little  or  no 
effect  when  the  predictions  are  averaged  over  an  interval  of  the 
range  variable. 

Consider  the  predictions  of  the  sound  field  by  both  ordinary 
and  modified  ray  theory  when  the  radiation  is  restricted  to  a ray 
bundle  dkv  centered  on  k^.  Both  theories  yield  similar  results, 
as  sketched  in  Fig.  4.  At  a fixed  receiver  depth  z , the  ray 
bundle  irradiates  a small  range  interval  on  the  way  down,  and 
another  small  interval  at  a slightly  greater  range  on  the  way  up. 
(These  correspond  to  the  terms  G*  and  G~  in  the  multipath  expan- 
sion for  the  monopole  field.)  Then,  at  a range  interval  equal 
to  the  period  X,  the  pattern  is  repeated.  (These  correspond  to 
the  terms  G++1  and  G~+1<) 

Differences  between  ordinary  and  modified  ray  theories 
are  exhibited  in  the  proportions  and  spacing  of  the  irradiated 
zones.  The  width  of  each  zone  is  proportional  to  the  range 
derivative,  being  given  by 

dr  = | 3r/8kv | dkv 

The  range  derivative  has  different  values  in  the  two  theories. 
Moreover,  the  squared  pressure  in  each  zone  is  inversely  pro- 
portional to  the  range  derivative  (Eqs.  18  and  37),  so  that  it 
also  differs  between  the  two  theories,  but  inversely  to  the  dif- 
ference in  zone  width.  The  result  suggests  that  we  should  con- 
sider the  range-averaged  prediction. 
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FIG.  4.  THE  PATTERN  OF  SOUARED  PRESSURE  AS  A FUNCTION  OF  RANGE,  WHEN  ONLY 
A SMALL  RAY  BUNDLE  IS  EXCITED.  THE  PATTERNS  PREDICTED  BY  ORDINARY 
AND  MODIFIED  RAY  THEORIES  ARE  SIMILAR  BUT  DIFFER  IN  DETAIL  (SEE  TEXT). 
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Following  procedures  earlier  found  to  be  convenient  in 
ordinary  ray  theory  [13],  we  shall  average  the  contribution  of 
each  ray  bundle  over  a range  Interval  equal  to  the  range  period  X, 
taking  all  ranges  to  be  equally  probable.  Then,  we  sum  these 
averaged  contributions  over  all  possible  ray  bundles.  For  sim- 
plicity of  exposition,  we  neglect  attenuation  and  set  |V  V^J  = 1. 

First,  we  need  an  expression  for  the  squared  pressure  ob- 
served in  the  zone  of  range  that  is  irradiated  by  a ray  bundle. 

By  the  same  algebra  that  was  used  earlier  for  monopole  radiation 
(Eqs.  33  through  35),  we  calculate  for  a general  source  the 
squared  pressure  observed  at  depth  z : 


p2  = pc  ( Zj  )I  = pcol/k  ( z j ) 


poj  3P/3kv  pco  3P/9kv 
k ( z ) 3A/3k  = k ( z ) 2nr|  3r/3kJ 

i V i V 

z 1 

where  I is  intensity,  P is  power,  and  dA  is  the  transverse  area 
of  a ray  bundle  dk^.  Now,  in  Sec.  5,  the  effective  source 
strength,  S of  any  source  was  defined  as  the  source  strength 

that  yielded  the  correct  value  for  power  in  a ray  bundle  when 
used  in  the  ordinary  ray-acoustical  formula  (Eq.  Al): 

3P/9kv  = 27rkvSeff/[pwkz  (z2 ) ] , 


where  z2  is  the  source  depth.  The  combination  of  these  results 
is : 


P2 (z,  ) = 


kv  Seff 


r kz(Zl)kz(z2)  |3r/3kvlz 


(5*0 
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To  form  the  average  contribution  of  any  ray  bundle  dk  , we 
multiply  p2  by  the  ratio  of  the  irradiated  zone  width  to  the 
range  period  X(kv),  a factor 

I 3r/3k  I dk 
v 1 v 

X 


The  result  is  seen  to  be  the  same  for  each  of  the  two  zones 
(corresponding  to  multipath  terms  G+  and  G~)  that  occur  in  each 
range  period.  (The  slight  variation  in  r is  ignorable  at  long 
ranges.)  Adding  these  two  contributions  and  summing  over  all 
ray  bundles,  one  finds  the  average  squared  pressure 


:P2 (z  j ) > = 


2 k 


v eff 


r X 


k (Zj )k  (z2  ) 


dk. 


(55) 


where  the  integral  extends  over  all  ky  such  that  the  modified 
ray  reaches  the  depth  z x (i.e.,  k2(z1)>0). 

It  is  often  convenient  to  introduce  as  a new  integration 
variable,  in  place  of  kv,  the  ray-acoustical  inclination  angle 
at  the  source,  02,  defined  by  Snell's  law: 

kv  = k ( z 2 ) cos02  , kz(z2)  = k(z2)  sin02  . (56  ) 


Then  Eq . 55  can  be  written  in  the  form 


<p2 (z  x ) > 


02>O 


2 S 

rr 


eff 
t an  8 ] 


de . 


(57) 


where  0j  is  the  inclination  angle  at  the  receiver's  depth  z2 
as  computed  by  Snell's  law.  This  integral  is  restricted  to 
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positive  02 , since  the  contributions  of  what  would  in  ordinary 
ray  acoustics  be  rays  at  ±02  are  both  included  in  the  effective 
source  strength.  The  contribution  of  k^  corresponding  to 
imaginary  values  of  02  can,  as  a practical  matter,  be  excluded 
without  significant  error  because  their  contribution  is  much 
weaker  than  that  of  the  real  angles. 

The  result  is  identical  in  form  with  the  corresponding 
result  of  ordinary  ray  theory  [Ref.  13,  Eq . 4].  However,  the 
value  of  the  integrand  differs  from  the  value  computed  from 
ordinary  ray  theory  for  two  reasons.  First,  the  range  period  X 
is  not  the  same  in  modified  as  in  ordinary  ray  theory.  Based 
on  the  illustrative  results  of  Murphy  and  Davis  [111  and  other 
data  presented  in  this  report,  this  error  appears  typically  to 
be  small  and  negligible.  The  value  of  X given  by  ordinary  ray 
acoustics  is  sufficiently  accurate  for  a range-averaged  modified 
ray  theory. 

The  second  difference  lies  in  the  value  of  the  effective 
source  strength  Se^.  As  shown  by  the  analytical  results  pre- 
sented elsewhere  in  this  report,  ordinary  ray  theory  can  lead 
to  serious  errors  in  evaluating  it. 

By  using  the  effective  source  strength,  found  by  a wave 
analysis  in  the  surface  layers  only,  methods  of  ordinary  ray 
theory  can  be  used  to  calculate  an  accurate  value  for  the  range- 
averaged  transmission,  although  the  point-to-point  transmission 
will  be  inaccurate.  Of  course,  there  is  no  need  to  perform  the 
point-to-point  calculations  when  the  averaging  can  be  done 
analytically,  with  calculations  based  on  Eq . 57- 
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7.  ISOGRADIENT  SURFACE  LAYERS 

In  this  section,  analytical  expressions  for  the  wave- 
theoretical  modifications  to  ray  theory  are  developed  for  the 
case  of  a source  in  an  isogradient  surface  layer.  It  is  assumed 
throughout  that  the  layer  is  deep  enough  that  ordinary  ray  (WKB) 
analysis  suffices  to  trace  the  spreading  of  sound  energy  below 
the  layer. 

Practitioners  of  wave  theory  will  have  guessed  that  the 
problem  we  shall  solve  is  that  for  which  the  Airy  functions  are 
precise  solutions  to  the  one-dimensional  wave  equation.  We 
postulate  an  environment  in  which  the  gradient  with  depth  of  k2 
is  a positive  constant: 

dk2(z)/dz  = B3  . (58) 

Since  k = m/c , this  prescription  implies  that  the  sound-speed 
gradient  is  negative  and  given  by 

g*  = -dc/dz  = iTfB3/k3  , (59) 

where  f = ui/2u  is  the  frequency.  Since  k varies  with  depth,  the 
sound-speed  gradient  is  not  strictly  constant  in  the  layer. 
However,  its  variations  are  negligible  in  practical  cases. 

We  require  solutions  to  the  one-dimensional,  separated 
wave  equation  (Eq.  8): 

d2g(z)/dz2  + k2 ( z ) g(z)  =0  , z > 0 , 

g(0)  = 0 

k2  = k2 (z)  -k2  , (60) 
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where  the  radial  wavenumber  k^  is  constant.  The  transformation 
of  coordinates 

y = B(z-zv) 

k * = B3(z-zv)  = B2y  , (61) 

brings  the  wave  equation  into  Airy's  form: 
d2g(y )/dy2  + y g(y)  =0  , y > y0 

g(y0)  = 0 > (62) 

where  y0  is  the  value  of  y corresponding  to  the  ocean  surface, 
z = 0.  Note  in  Eq . 6l  that  z is  the  depth  corresponding  to  the 
ray  vertex,  or  turning  point,  where  kz  vanishes.  Its  value  may 
be  positive  or  negative,  corresponding  in  the  latter  case  to  the 
depth  at  which  the  ray  would  have  vertexed  if  the  layer  had  ex- 
tended above  the  actual  sea  surface. 

The  solution  to  this  boundary  valve  problem  is  well  known: 

g(y)  = C M(y ) sin [ ^ y ) ] 

<P0  = >My0)  , (63) 

where  C is  an  arbitrary  constant  and  M and  \p  are  real  functions 
defined  In  terms  of  Airy  functions: 

M'y)ei'My'>  = Ai(-y)  + i Bi(-y)  . (64) 
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Here,  Ai  and  Bi  are  the  standard  Airy  functions,  for  which  various 
analytical  properties  and  tables  of  values  are  found  in  standard 
texts  [Abramowitz  and  Stegun,  Ref.  14] . Equation  63  is  the  linear 
combination  of  Airy  functions  that  matches  the  boundary  condition 
at  y = y0. 


Normalization 


In  the  preceding  development  of  modified  ray  theory,  the 
solution  to  the  one-dimensional  wave  equation  was  normalized 
in  a very  particular  way.  Specifically,  the  desired  solution 
is  the  response  to  an  upwardly  traveling  wave  that  has  unit 
amplitude  and  zero  phase  at  a reference  depth  z . (See  Eqs.  7 
and  9.)  The  reference  depth  is  chosen  to  be  sufficiently  deep 
that  the  asymptotically  precise,  ordinary  WKB  expression  ac- 
curately describes  the  local  field,  but  it  is  otherwise  arbitrary. 
That  is,  the  constant  C in  Eq.  63  is  to  be  chosen  so  that,  for 
z and  Zj  both  deep  in  the  asymptotic  region, 

2 kz  ^ z 1 ^ 

g = g ( z , z ! ) = n * 

U k2(z) 


1 

•z 

4 <J>  + 

k dz 

2 ^u 

z 

l z,  J 

(65) 


(Compare  Eq . 21.)  Here,  <J>u  has  the  significance  that 

1 4> 

Vu(z1 ) = -e  u (66) 


is  the  plane-wave  reflection  coefficient  of  the  ocean,  0<z<z1} 
as  observed  at  the  depth  z, . 
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The  asymptotic  forms  of  M and  ip  for  large  positive  y are 
[Ref.  14;  Eqs . 10.4.78,  10.4.79]: 

-k  _k 

M ~ it  2 y 14 

>4  ~ ^-TT  - (2/3)y3^ 

dip/ dy  ~ -y'2  , (67) 

which  are  accurate  for  y3  >>  5/32.  These  forms  can  be  translated 
by  means  of  Eq . 6l  into: 

M ~ (B/irk  )'2 
• z 

My)  ~ 4;  k dz 

^ x = My  j ) , (68) 


where  y t is  the  value  of  y corresponding  to  z = Zj.  Now,  a direct 
comparison  of  the  desired  form,  Eq . 65,  with  the  general  Airy 
function  solution,  Eq . 63,  yields  the  required  value  of  the  con- 
stant : 


4>u  = arg[-Vu(z1)]  = 2(ipg-ip1) 


(69) 
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Error  in  Range  Period 

In  this  subsection,  we  evaluate  the  error  in  range  period  X 
that  arises  from  the  use  of  ordinary  ray  acoustics. 

The  range  period  is  related  to  the  phases  of  the  plane-wave 
reflection  coefficients  for  the  upper  and  lower  halves  of  the  ocean 
(Eq.  32).  In  this  study,  we  are  concerned  with  the  error  arising 
from  the  upper  half.  The  error  in  the  range  period  Is  due  to 
the  difference  between  the  wave-theoretical  phase  and  the  asymp- 
totic expression  for  that  phase,  since  the  asymptotic  expression 
is  that  used  in  ordinary  WKB  analysis.  The  phase  of  the  re- 
flection coefficient  in  the  present  isogradient  problem  is  given 
in  Eq . 69-  Only  the  term  \pQ  = ip(yg)  differs  from  its  asymptotic 
form,  since  the  observation  depth  is  by  definition  so  large 
that  ip  j is  given  by  its  asymptotic  value.  Thus,  we  are  led,  by 
the  combination  of  Eqs.  32  and  69,  to  this  expression  for  the 
error  in  the  range  period  (i.e.,  the  ray-acoustical  value  less 
the  wave-acoustical  value): 

AX  = -2  3 [ ipas  ( y 0 ) - iMy0)]/3kv  , (70) 

O Q 

where  tp  denotes  the  asymptotic  form  of  ip , and  y0  is  implicitly 
a function  of  the  ray  parameter  ky. 

Two  asymptotic  forms  must  be  distinguished,  depending  on 
the  sign  of  y 0 : 

1 tt/4  - 2y  3/2/3  for  y > 0 
<P' (y)  = < 

) tt/2  for  y < 0 . (71) 
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The  first  has  been  given  in  Eq . 67.  The  second  follows  from  the 
known  asymptotic  behavior  [74]: 

Ai(x)/Bi(x)  ~ 0 as  x ~ 00  . 


The  derivative  in  Eq.  70  can  usefully  be  transformed  by  successive 
application  of  Eqs.  60,  6l,  and  59: 


? k 

3 „ , 3 „ , 3 _ * v d 

tt-, — = 2 k = -2  k,  = — 

9kv  V 3(kJ)  3(kJ)  B2  dy 

Then,  the  error  in  range  period  takes  the  form 

y)  . y,  > °|  | 

0 . .vD  < 0 1 


AX 


A k 


B" 


I ib  ' 


^ ' (y  0 ) + 


(72) 


where  ip  ' = dip/dy . 

The  environmental  variables  and  ray  parameters  are  rather 
intricately  interwoven  in  the  variables  and  parameters  that  appear 
in  Eq.  72.  However,  it  is  possible  to  show  that  the  error  in 
wavelengths  is  a function  of  only  two  variables:  (i)  the  ratio 

of  frequency  to  sound-speed  gradient,  and  (ii)  a ray-theory  angle 
at  the  surface  defined  by 

kv  = k(0)  cos60  . (73) 


This  angle  is  imaginary  for  rays  that  vertex 
although  it  is  still  a useful  variable,  being 
depth  of  the  turning  point  expressed  in  units 
Graphs  of  typical  values  are  shown  in  Figs.  5 


below  the  surface 
then  related  to  the 
of  the  wavelength, 
and  6 . 


h6 


SURFACE  REFLECTED  RAYS 


FIG.  6.  ERROR  6X  IN  THE  RAY-THEORETICAL  RANGE  PERIOD  AS  A FUNCTION  OF  Isinel  AT  THE  SURFACE, 
FOR  (f/q*)  = 50  THE  LEFT  HALF  OF  THE  GRAPH  APPLIES  TO  ORDINARY  RAYS  THAT  VERTEX  AT 
A DEPTH  z BELOW  THE  SURFACE,  FOR  WHICH  0 IS  IMAGINARY. 


Report  No.  3326 


Bolt  Beranek  and  Newman  Inc. 


It  will  be  noted  that  the  largest  error  in  the  range  period 
occurs  for  the  ray  whose  turning  point  is  at  the  surface.  This 
maximum  error  is  given  by 


AX 


max 


1.261 

fl/3  r*2/3 


km 


(74) 


where  f is  in  Hz,  g*  in  sec  1 , and  the  value  1462  m/s  has  been 
used  for  sound  speed.  This  relation  is  plotted  in  Fig.  7. 


Equation  74  and  Fig.  7 are  both  very  deceptive.  Two  cautions 
must  be  stressed.  First,  the  analysis  assumes  that  the  constant- 
gradient  layer  is  infinitely  deep.  It  is  applicable  to  a layer 
of  finite  depth  only  if  the  ray  in  question  (horizontal  at  the 
surface)  steepens  within  the  layer  to  an  extent  that  asymptotic 
analysis  is  valid  at  greater  depths.  Taking  values  of  y > 1 as 
defining  the  asymptotic,  ordinary  WKB  region,  one  can  calculate 
the  layer  depth  required  for  validity  of  the  results.  The  limiting 
contours  in  the  case  of  layer  depths  of  50  and  100  m have  been 
plotted  in  Fig.  7. 


The  second  caution  concerns  the  change  in  range-period  error 
with  changes  in  the  gradient,  frequency  being  held  constant. 

Figure  7 shows  the  error  increasing  as  the  gradient  decreases ! 

The  result  seems  to  contradict  our  intuition,  which  is  based  on 
the  expectation  of  small  errors  in  ray  theory  when  the  gradient 
is  small.  The  paradox  is  resolvable.  First,  the  limit  as  gradient 
approaches  zero  is  not  pertinent  to  physical  problems,  because  of 
the  limited  depth  of  actual  isogradient  layers.  Second,  the  rate 
of  change  of  error  with  gradient  depends  on  the  steepness  of  the 
ray,  as  a comparison  of  Figs.  5 and  6 will  indicate.  Only  for  a 
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50 


FIG.  7.  CONTOURS  OF  MAXIMUM  ERROR  IN  THE  RAY-THEORETICAL  RANGE  PERIOD.  THE  MAXIMUM  OCCURS 
FOR  THE  RAY  HORIZONTAL  AT  THE  SURFACE.  THE  RESULTS,  VALID  FOR  VERY  DEEP  LAYERS,  ARE 
VALID  ONLY  TO  THE  RIGHT  OF  THE  LIMITS  SHOWN  FOR  FINITE  LAYER  DEPTHS. 
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narrow  beam  of  rays,  centered  on  the  ray  that  grazes  the  surface, 
does  the  error  grow  with  decreasing  gradient.  The  width  of  that 
beam  becomes  smaller  as  gradient  decreases  or  as  frequency  in- 
creases. Eventually,  the  beam  width  becomes  insignificant. 
However,  Figs.  5 and  6 demonstrate  that  a significant  beam  has 
this  anomalous  behavior  at  low  frequencies. 


Effective  Source  Strength  of  Monopole 

The  general  formula  for  effective  source  strength  of  a 
monopole  (volume  source)  is  given  in  Eq . 44.  For  a unit  monopole 
(SQ  = 1),  the  effective  source  strength  D0  involves  the  solution 
to  the  separated  wave  equation,  g^ , evaluated  at  source  depth  z2. 
The  suitably  normalized  function  g for  the  isogradient  environment 
is  given  by  Eq . 63  with  the  constant  C given  in  Eq.  69.  The 
combination  of  these  results  can  be  written  as 

D „ = 4ir[y^M2  sin2  ( ^ Q ) ] , (75) 

y~y  2 


where  y2  is  the  value  of  y at  z = z2.  (The  explicit  appearance 
of  the  environmental  parameter  B has  been  eliminated  from  Eq.  75 
by  use  of  Eq.  6l.) 


The  variable  y is  dependent  on  frequency,  sound-speed  gradient, 
and  a variable  identifying  the  ray.  It  is  convenient  to  use  as 
ray  variable  the  inclination  angle  02  at  source  depth,  as  given 
by  Snell's  law  (Eq.  56).  Then,  using  Eqs.  59  and  6l,  one  can 
show  that 

y 2 = ^(z2)/B2  = (TTf/g*)273  sin20  2 , (76) 


where  g*  is  the  negative  of  sound-speed  gradient  and  f is  frequency. 


51 


Report  No.  3326 


Bolt  Beranek  and  Newman  Inc. 


This  wave-theoretical  value  of  D0  is  to  be  compared  with  the 
ordinary  WKB,  ray-theoretical  value  compared  from  Eq . 45,  and  with 
Lloyd's  mirror  formula,  Eq . 46.  Figure  8 shows  this  comparison 
for  typical  values  of  the  parameters.  It  is  interesting  to  see 
that  a precise  application  of  ordinary  ray  theory  is  less  accurate 
than  Lloyd's  mirror  formula.  The  wave-theoretical  result  (labeled 
"modified  ray  theory")  indicates  that  the  effective  source  strength 
is  significantly  larger  than  the  approximations  for  it,  when  the 
angle  02  is  small. 

The  general  functional  dependence  of  DQ  on  02,  involving 
the  functions  M(y)  and  ip(y)  , is  very  intricate.  However,  a simple 
approximation  can  be  derived  for  the  cases  of  principal  interest. 
First,  we  are  concerned  with  sources  close  to  the  ocean  surface, 
in  the  "decoupling  layer",  where  the  sine  function  in  Eq . 75  is 
small  and  its  argument  is  small  relative  to  tt/2.  Indeed,  one  can 
show  that  the  asymptotic,  ordinary  WKB  solution  is  good  when  this 
is  not  the  case.  Then,  the  sine  is  approximately  equal  to  its 
argument,  which  can  be  rewritten  using  the  law  of  the  mean: 

sin(t|>o-t|j)  * t|/0  - tp  = i|/'(y)  (yQ-y) 

where  \p ' is  the  derivative  of  ip  and  y is  a mean  value  of  the 
argument,  y(  < y < y . But  the  derivative  of  the  phase,  ip  ’ , is 
related  to  the  amplitude  M by  an  identity  for  Airy  functions 
[Ref.  14,  Eq.  10.4.71]: 

*p'  (y)  = -1/nM2  ( y ) 

We  introduce  these  approximations  into  Eq . 75*  Further,  (y0-y2) 
is  related  to  zz  by  Eq . 6l,  and  y"2  to  sin02  by  Eq . 76.  Finally, 
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FIG.  8 


. EFFECTIVE  SOURCE  STRENGTH  D OF  A UNIT  MONOPOLE  SOURCE,  CALCULATED 
FROM  EQ.  75.  0 


Report  No.  3326 


Bolt  Beranek  and  Newman  Inc. 


we  introduce  the  approximation 


M2 ( y )/M4 (y ) = 1/M2(0)  = 1.983 


which  is  accurate  when  both  y2  and  yQ  are  small  - that  is,  in  the 
region  where  the  wave-theoretical  correction  is  important.  The 
resulting  approximation  for  effective  source  strength  is 
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/M! 

\t!  \t ) 


sin6 . 


(77) 


The  range  of  validity  of  the  approximation  is  derived  below  (Eq.  79) 

This  approximation  is  to  be  contrasted  with  a similar  approxi- 
mation to  Lloyd's  mirror  formula  (Eq.  46): 


= 4 sin2(k  z ) « 4 ( k z )2 

0 Z 2 Z 2 


or,  writing  kz  = k sin©2. 


= 158(z2/A)2  sin202 


(78) 


The  result  is  valid  when  D 


LM 


• 0 <<  4,  which  is  not  a serious  limi- 

tation in  the  cases  of  present  interest. 


The  striking  difference  between  Eqs.  77  and  78  lies  in  their 
dependence  upon  02 . This  difference  causes  the  error  in  the 
Lloyd's  mirror  formula  to  grow  with  decreasing  values  of  inclina- 
tion angle  02.  Figure  9 shows  the  results  of  precise  calculations 
of  the  error  as  a function  of  02  and  the  ratio  of  frequency  to 
gradient.  Also  plotted  on  the  figure  is  the  estimate  of  the 
error  formed  by  the  ratio  of  Eq . 78  and  Eq . 77.  Within  the  pre- 
cision of  these  approximations,  the  error  is  independent  of  source 
depth . 
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FIG.  9.  ERROR  IN  THE  EFFECTIVE  SOURCE  STRENGTH  OF  A UNIT  MONOPOLE  SOURCE  THAT  RESULTS 

FROM  LLOYD'S  MIRROR  FORMULATION.  THE  DASHED  LINES  ARE  APPROXIMATIONS  CALCULATED 
FROM  EOS.  77  AND  78. 
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Finally,  the  limit  of  validity  of  the  approximate  wave- 
theoretical  result,  Eq.  76,  can  be  established  from  the  observa- 
tion that  the  approximation  is  good  when  it  exceeds  the  Lloyd's 
mirror  formula.  Hence,  Eq . 77  is  valid  if 

gVf  > 12.5  s i n 3 0 2 , ( 79  ) 

which  is  equivalent  to  y2  < 0.40  (see  Eq . 77).  When  this  in- 
equality fails,  the  Lloyd's  mirror  formula  is  accurate,  as  is 
shown  by  the  convergence  of  the  precise  values  in  Fig.  8.  When 
this  inequality  holds,  both  the  Lloyd's  mirror  formula  and  the 
precise  result  of  ordinary  ray  theory  are  inaccurate. 

Effective  Source  Strength  of  Vertical  Dipole 

The  general  formula  for  the  effective  source  strength  of  a 
vertical  dipole  (force  source)  is  given  in  Eq.  51.  For  a unit 
force  F,  the  effective  source  strength  involves  the  spatial 

derivative  of  the  solution  to  the  separated  wave  equation, 
evaluated  at  source  depth.  With  the  introduction  of  the  solution 
for  the  isogradient  environment  (Eqs.  63  and  69)  into  Eq . 51,  the 
effective  source  strength  per  unit  force  can  be  written  as 

DVF  = 47T^kz^Z2  )/k(z2  [y~^N2sin2  ( ^ „ -h  ) ]y_y  , (80) 

where 

N(y)ein(y)  = Ai’(-y)  + i Bi'(-y)  . (8l) 

This  wave-theoretical  value  of  D^p  is  compared  with  the 
ordinary  WKB,  ray-theoretical  value  (Eq.  52)  and  the  Lloyd's 
mirror  formula  (Eq.  53)  in  Fig.  10.  As  was  the  case  for  the 
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DEPRESSION  ANGLE  AT  SOURCE  (DEG) 


FIG.  10.  EFFECTIVE  SOURCE  STRENGTH  Du(-  OR  A VERTICAL  DIPOLE,  CALCULATED 
FROM  EQ.  80. 
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monopole  source,  the  wave-theoretical  value  of  source  strength 
(labelled  "modified  ray  theory")  is  significantly  larger  than  the 
approximations  when  the  ray  inclination  angle  is  small. 

An  approximation  to  Eq . 80  displays  the  explicit  dependence 
on  physical  parameters.  One  starts  with  the  expansion 

<1*0  - n = - (n-ijO 

Numerical  experience  indicates  that  the  first  term  on  the  right 
is  small  and  negligible  with  respect  to  the  second  term  in  the 
range  of  parameters  of  interest.  The  second  term  ranges  from 
tt/3  to  tt/2  for  all  y > 0 (02>O).  But  the  first  term,  which 
appears  in  Eq . 75  for-  the  monopole,  is  small  relative  to  tt/2  in 
the  surface  decoupling  layer.  The  sine  of  the  second  term  is 
now  transformed  by  an  identity  for  Airy  functions  [Ref.  14,  Eq . 
10.4.74], 


sin  (h-ifO  = 1/tt  MN  , 
with  the  result 

DyF  = 4[kz(z2 )/k(z2  ) ]2  [n  y2M2(y2)]_1 

= 4 sin2  02  [iry2M2(y2)]-1  . (82) 

It  is  interesting  to  note  the  asymptotic  behavior  of  Eq.  82. 
From  the  asymptotic  expression  for  M (Eq.  67),  which  is  a good 
approximation  for  yz  1,  one  finds 
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DyF  s 4 sin2 0 2 

which  is  identical  with  an  approximation  of  the  Lloyd's  mirror 
formula  (Eq . 53 ) : 

= 4 sin202  cos2(kzz2)  = 4 sin202  . (83) 

(The  approximation,  good  for  kzz  small,  is  not  an  important 

limitation  in  the  practical  problems  of  interest.)  Hence,  we 

LM 

have  demonstrated  analytically  the  confluence  of  Dyp  and  Dyp  for 
y2  > 1,  which  is  observed  in  the  numerical  data  of  Fig.  10. 

The  asymptotic  form  for  M2  becomes  significantly  poor  (errors 
of  1 dB  or  more)  for  values  of  y2  < 0.5.  However,  the  value  of  M2 
does  not  vary  widely  in  the  interval  0.5  > y2  >_  0,  so  that  the 
geometric  mean  value 

M2  = [M2(0)  M2 (0 . 5 ) = 0.4  32 

can  be  used  without  introducing  errors  greater  than  0.7  dB.  Using 
this  approximation  in  Eq . 82  and  introducing  the  physical  para- 
meters, we  find  the  approximation 

Dyp  s 2.01  (g*/f)1/*  sin02  . (84) 

This  approximation  is  good  for  real  02  so  long  as  it  exceeds  the 
asymptotic  expression  previously  shown  to  be  the  same  as  Eq . 83. 
Hence,  Eq . 84  is  applicable  if 

(gVf)  > 7.9  sin 3 0 2 , ( 85  ) 

which  is  equivalent  to  y2  < 0.54  (see  Eq.  77). 
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As  in  the  monopole  case,  we  find  from  wave  theory  that  the 
effective  source  strength  varies  as  sin©  when  0 approaches  zero. 
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